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We show that the energy spectrum of the Bose-Hubbard model amended by a static field exhibits 
Wigner-Dyson level statistics. In itself a characteristic signature of quantum chaos, this induces the 
irreversible decay of Bloch oscillations of cold, interacting atoms loaded into an optical lattice, and 
provides a Hamiltonian model for interaction induced decoherence. 
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The Bose-Hubbard Hamiltonian serves as a paradigm 
in the field of quantum phase transitions pj. Recently, 
this model was realized in experiments on ultracold 
atoms loaded into a three-dimensional optical lattice Q , 
opening new perspectives for the laboratory study of cor- 
related bosonic systems. Consequently, new theoreti- 
cal work on the Bose-Hubbard model was stimulated, 
which, in particular, addresses the response to a static 
field 0, 0, 01 - a question which shifts the focus from the 
Bose-Hubbard ground state (which is mostly studied in 
the literature) to dynamical and spectral properties of the 
system. In single-particle quantum mechanics, these are 
associated with Bloch oscillations, in the time domain, 
and with the emergence of a Wannier-Stark ladder, in 
the energy domain 

The present Letter is devoted to the spectral properties 
of the Bose-Hubbard Hamiltonian under the additional 
action of a static field or, equivalently, to the Wannier- 
Stark problem for interacting bosons. Our analysis is 
formulated in a spirit close to ongoing experiments on 
cold atoms in optical lattices 0, > and we assume that 
the atoms are in the "super-fluid phase" , i.e. they are de- 
localized over the lattice in the absence of any external 
perturbation. This latter assumption distinguishes the 
present work from previous contributions 0, 3 devoted 
to the Mott insulator phase, and restricts the values of 
the hopping matrix element J and of the on-site inter- 
action energy W to the range W/J < 5.8 (see and 
references therein). To be specific, we fix J = 0.038 and 
W = 0.032 - the experimental values (in units of pho- 
ton recoil energies) for rubidium atoms in optical lattices 
with a potential well depth of approx. ten photon recoils 
- throughout the sequel of this Letter. With a], hi, and 
iii the single particle creation, the single particle anihila- 
tion and the number operator at site / of the lattice, the 
total Hamiltonian reads: 
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FIG. 1: Energy spectrum (top left) together with Wannier 
state projections |c;| 2 of the field-free k = eigenstate of Q, 
for increasing values of the static field < F < 0.1. N — 1 
particle is loaded into a lattice with eleven sites (L = 5). The 
level dynamics illustrate the transition from a Bloch-band to 
a Wannier ladder, associated with the progressive localisation 
of the wave function in I. 



where the strength of the static field F or, more precisely, 
the Stark energy (the period of lattice is set to unity) will 
be our free parameter. 

Let us first address the issue of boundary conditions. 
It is well known that, for an infinite lattice, there is no 
smooth transition between the spectrum at F = and 
F =^ 0. Formally, this is due to the fact that for any non- 
vanishing value of F the Hamiltonian Q is an unbounded 
operator, whereas it is bounded for F = 0. However, for 
a lattice of finite size, —L < I < L, the operator is 
always bounded and, hence, the spectrum of the system 
changes continuously as a function of F, as illustrated 
by the numerically generated level dynamics in the top 
left panel of Fig.^ for N = 1 and Dirichlet (i.e., vanish- 
ing) boundary conditions, to L — 1). As F is increased, 
the spectrum evolves from a Bloch spectrum with ener- 
gies E(k) = — J cos[-Kk/{2L + 2)], k — —L, . . . , L, into a 
Wannier-Stark ladder Ei ~ Fl, I = —L, . . . , L. The other 
panels in Fig. ^ show the evolution of the field- free k = 
eigenstate in the basis of the Wannier states, with in- 
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creasing F |8j . The progressive localisation of the atomic 
wave function in I with F is known as Stark localisation. 

When discussing the time evolution of a wave function 
governed by 0J, it is preferable to use periodic bound- 
ary conditions instead of Dirichlet. To do so, one first 
eliminates the static term in by transforming to the 
interaction representation, where the hopping and the 
on-site term in Q define the unperturbed Hamiltonian, 
hence 
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and then identifies the site I — L + 1 of the lattice with 
I = —L. This choice has the advantage that the time 
evolution operator of a system of non-interacting atoms 
over one Bloch period Tb — 2nh/F coincides with the 
unit matrix, independently of the size of the system. This 
facilitates the analysis of the dynamics in the thermody- 
namic limit N, L — * oo J5J- In what follows, we shall use 
Dirichlet boundary conditions when calculating eigenval- 
ues, and periodic boundary conditions when simulating 
the dynamics 

Our analysis of the spectrum of the multi-particle sys- 
tem follows the one for the single-particle problem. 
Let us assume for the moment that there are no atom- 
atom interactions, i.e. W = 0. As already mentioned, 
for large values of F the single-atom energy levels form 
a Wannier ladder, and the energies of an iV-atom system 
are consequently given by 
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where the m; (m = to-l, . . . , tul, 1111 = -^0 are ^ ne 
occupation numbers of the Wannier-Stark states. Note 
that, in general, many different sets m correspond to the 
same total energy, and the TV-particle Wannier ladder- 
levels E m = Fltot are, thus, typically degenerate. The 
./V-particle wave function associated with a given level 
E m can be constructed from single-particle Wannier- 
Stark states \ipi) by an appropriate symmetrisation pro- 
cedure. In the basis of Fock states (symmetrised products 
of Wannier functions \n}), an arbitrary Wannier-Stark 
state, at finite F, is given by the sum 
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and in the limit F — > oo only one coefficient Cn"^ with 
n = m differs from zero in Eq. Q. On the contrary, 
in the opposite limit F — > 0, almost all expansion co- 
efficients are non-zero and the Wannier-Stark states ap- 
proach TV-particle Bloch states with (once again, degen- 




FIG. 2: Energy spectrum of the Hamiltonian Q as a function 
of the static field F, for particle number N = 3 and lattice 
size 2L + 1 = 11. Particle interaction strength W = 0.032, 
hopping matrix element J = 0.038. The inset zooms into the 
central part of the spectrum, in the range < F < 0.01. 
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the straightforward TV-particle generalization of the 
above one-particle result. 

Let us now include the effect of atom-atom interac- 
tions. Figure |21 shows the energy levels of the Hamil- 
tonian JTJ as a function of F, for N — 3 atoms loaded 
into a lattice with 11 sites (i.e., L = 5). As expected, 
the atom-atom interactions remove the above-mentioned 
degeneracy - for small F the spectrum appears dense (al- 
most continuous), and for large F the degenerate levels 
of the Wannier ladder split into "Wannier-ladder energy 
bands" (see Eq. JHJ) below). In this latter limit, the spec- 
trum and the associated Wannier-Stark states can still 
be found analytically. Indeed, since the hopping term in 
Eq. couples only those Fock states separated by one 
single quantum in the Stark excitation, one has 
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The perturbative results I|tjl7|) cannot hold when F < 
J. Moreover, the complex level dynamics which are 
borne out for small F in the inset in Fig. [21 indicate 
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FIG. 3: Solid line: Cummulative nearest neighbour level spac- 
ing distribution J(s) = J P(s')ds' for normalised spacings 

s = AE/AE (AE is the average level spacing in the central 
part of the spectrum), with F — 0.01 (top) and F = 0.04 (bot- 
tom). N = 4 atoms loaded into a lattice of size 2L + 1 = 11. 
The dashed and dash-dotted lines indicate GOE and Poisson 
cummulative distributions, respectively. 



that any attempt to assign a set of quantum numbers 
to individual levels is bound to fail for F < J. In- 
stead, a statistical analysis of the spectrum is appropri- 
ate in this situation. For that purpose, the upper part of 
Fig. 13 presents the cummulative distribution of the spac- 
ings between adjacent energy levels, for N — 4 atoms 
loaded into a lattice with eleven sites, at F = 0.01 [To|. 
Clearly, the normalised energy intervals s = AE/AE ex- 
hibit GOE (Gaussian Orthogonal Ensemble) statistics, 
P{s) = (ir 2 /6)s exp(— irs 2 /4), a hallmark of quantum 
chaos Thus, for weak static fields, the system (Q can 
be regarded as a quantum chaotic system. The origin of 
"quantum chaos", i.e., of the strongly F-dependent, non- 
perturbative mixing of energy levels can be understood 
here as a consequence of the interaction induced lifting of 
the degeneracy of the multi-particle Wannier-Stark levels 
in the cross-over regime from Bloch to Wannier spectra, 
making nearby levels strongly interact, for comparable 
magnitudes of hopping matrix elements and Stark shifts. 
In contrast, for large F (and in the limit of large L,N), 
the nearest neighbour distribution tends towards Pois- 
sonian statistics, P(s) = exp(— s), as evident from the 
lower part of Fig. [3] 

Which are the physical manifestations of the irreg- 
ular spectrum of QJ, at intermediate field strengths? 
To answer this question, we consider the Bloch oscil- 
lations of the mean atomic momentum, which can be 
observed rather easily in state of the art experiments 
0, ^3- I n the absence of atom-atom interactions, the 
average momentum p(t) of the atoms oscillates with the 
Bloch frequency lob — F/Ti. As shown in 5], the pres- 



ence of atom-atom interactions modifies the Bloch dy- 
namics, and p(t) exhibits an additional beating signal at 
frequency low — W/%: 

p(t) = tr \pp(t)] = NJf(t) wa.(u B t) , (8) 
f(t) = exp (-2n[l - cos(wv^)]) , 

with n the filling factor, i.e., the average number of atoms 
per lattice cite, and pit) the single particle density matrix 
with elements (in the Wannier state basis) 

Pl>v {t) = <tt(t)$6H*(t)> ■ (9) 

The appearance of the new frequency low originates in 
the splitting of the Wannier ladder levels into "energy 
bands" - see our above discussion. It must be stressed 
that the result JSJl is valid only for large values of the 
static field, where the spectrum is regular. Consequently, 
it is to be expected that for weak static fields the atomic 
Bloch oscillations will be qualitatively different, due to 
the irregular/chaotic structure of the spectrum. Indeed, 
numerical simulations of the dynamics indicate that, in 
the weak field regime, the Bloch oscillations decay ir- 
reversibly on rather short time scales. As an example, 
Fig. 0] shows the behaviour of the scaled (p — > p/NJ) 
momentum for F = 0.05, N = 7, L = 3 (top), and 
N — 9, L = 4 (bottom) 13]. After only few Bloch peri- 
ods, the mean momentum has almost comletely decayed 
but for feeble residual fluctuations, which rapidly van- 
ish as the Hilbert space dimension is increased (compare 
top to bottom). Thus, for a weak static field, the enve- 
lope function in Eq. © approaches f(t) ~ exp(— tjr) in 
the thermodynamic limit, which appears to be reached 
rather quickly with increasing system size, by virtue of 
the results shown in Fig. Note that the decay of the 
Bloch oscillations is due to the decay of the off-diagonal 
elements of the one-particle density matrix J^J, and that, 
hence, the time r can equally be considered as the de- 
coherence time for a system of interacting bosons. The 
dependence of r on the system parameters hitherto re- 
mains an open problem. 

Finally, let us briefly discuss the conditions for the 
observed chaos-transition in the Bose-Hubbard model. 
Our numerical simulations of the system dynamics, per- 
formed for fixed ratio W/ J and different values of N and 
L (0.2 < n < 1.2, L < 10, N < 10), suggest the condition 

SI ~ n" 1 , (10) 

as a criterion of the transition to chaos, where SI denotes 
the localisation length of the single-particle wave func- 
tion on the lattice (51 ~ J/F for F < J, and SI ~ 1 
for F > J), and n" 1 as the inverse filling factor has 
the meaning of an average particle distance. It is clear, 
however, that condition (|10fl cannot be universal, since 
it does not account for the on-site energy W . Indeed, 
for W — > 0, the particle-particle interaction vanishes, 
and the system is integrable for arbitrary F. On the 
other hand, when W oo, the Bose-Hubbard model 
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FIG. 4: Irreversible decay of atomic Bloch oscillations in the 
presence of a weak static field F = 0.05, for a filling factor 
n = 1 and lattice sizes 2L + 1 = 7 (top) and 2L + 1 — 9 
(bottom). Comparison of both plots suggests that the resid- 
ual fluctuations of p(t) rapidly die out as the system size is 
increased. (The respective dimensions of Hilbert space are 
N = 1716 (top) and TV = 24310 (bottom).) 



can undergo a Mott transition into the insulating phase, 
where its response to the static field has a very differ- 
ent (resonant-like) character 09- It therefore remains 
a challenging theoretical problem to formulate general 
criteria for the chaos-transition. 

To conclude, we have shown that the spectrum of the 
Bose-Hubbard Hamiltonian amended by a static field 
(and at fixed particle-particle interaction corresponding 
to the "superfluid" regime in field-free case) is cither 
regular or irregular, depending on the relative strength 
of the hopping matrix element and the external pertur- 
bation. In particular, wc have seen that the irregular 
level structure at intermediate strengths of the static field 
manifests in a rapid decay of the Bloch oscillations of the 
mean atomic momentum, and that the time scale of this 
decay provides a direct measure for the decay of particle- 
particle coherences accross the lattice. Hence, chaotic 
dynamics of cold, interacting atoms loaded into a one di- 
mensional optical lattice allow for experimental probing 
and control of interaction induced decoherence. 

We thank Boris Fine and Hcnning Schomerus for useful 
discussions. 
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